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We study the entanglement dynamics of two non-interacting, spatially 
separated qubits subject to local environment noises. Based on exactly 
solvable models for non-Markovian amplitude damping and phase damp- 
ing noises, we are able to analyze the entanglement dynamics of the two 
qubits for different coupling bandwidths and different detunings. We show 
that entanglement oscillations can occur for both amplitude and phase 
damping noises. Moreover, we demonstrate that changing the coupling 
bandwidth can lead to crossover between dissipative and non-dissipative 
entanglement dynamics, while varying the detuning controls the crossover 
between strong and weak coupling limits. Our findings can help provide a 
synthesized picture for the entanglement dynamics of two qubits subject 
to local environment noises. 
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I. Introduction 

Advancements of quantum information sciences during the past few decades 
have been bringing new lights into researches on the fundamental issues of 
quantum mechanics pQ. Among them, studies of entanglement dynamics for 
open quantum systems have occupied a central part due to the crucial role of 
quantum entanglement in many aspects of quantum information sciences [2]. 
For general multi-qubit systems, the study of entanglement dynamics remains a 
major challenge due to the lack of a manageable quantitative measure for entan- 
glement in such systems [3]. In the case of two-qubit systems, however, there 
have been reliable and easily computable entanglement measures which have 
greatly assisted investigations of entanglement dynamics in these systems. In 
particular, the degradation of entanglement in two-qubit systems in the presence 
of environment noises has been an issue of great theoretical and experimental 
interests [I]. 

In this paper, we shall examine the entanglement dynamics of two spatially 
separated, non-interacting qubits subject to local environment noises. This may 
correspond to, for instance, a pair of entangled two-level atoms located in two 
different (lossy) cavities, where each atom is separately coupled to the cavity 
modes. In order to study the entanglement evolution of the qubits under the 
action of environment noises, we will consider the cases of amplitude damping 
noise, which can cause energy relaxations in the qubits, and of phase damping 



noise, which can degrade the phase coherence of the qubits pQ. The Marko- 
vian dynamics of the qubit entanglement for such systems have been studied 
extensively [4] . Novel effects such as the complete destruction of entanglement 
within finite time (the entanglement sudden death) have been predicted and 
observed experimentally [U [5J [7] . There have also been researches on the 
non-Markovian entanglement dynamics in such systems. In particular, entan- 
glement oscillations due to the finite memory time of the environment coupling 
have been predicted theoretically for amplitude damping noise [9] , and ob- 
served experimentally in photonic systems [ID]. In this paper, based on analyti- 
cally solvable models that exhibits non-Markovian characteristics, we will show 
that for both amplitude damping and phase damping noises, crossovers between 
oscillatory and monotonic decaying entanglement evolution can occur in con- 
nection with the transition between the non-dissipative and dissipative limits of 
the environment coupling, which may also be regarded as a crossover between 
the non-Markovian and the Markovian limits. Furthermore, our model shall 
also enable us to examine the crossover between the limits of weak and strong 
environment couplings. Therefore, our results demonstrate not only that entan- 
glement oscillations can occur also for two qubits under phase damping noises 
(cf. Ref. [H]), but also offer detailed analysis for the entanglement dynamics 
of two qubits in various regimes which can help provide a synthesized view for 
the related problems. It should be noted, however, that here our emphasis will 
lie in the entanglement dynamics of the qubits, rather than the different mea- 
sure^) for the "non-Markovianity" of the dynamical maps that we will discuss 
[El H31 E] ; even though there are connections between the two (see later). In 
other words, it is the physical reasons behind the appearance/disappearance of 
entanglement oscillations that will concern us, instead of any specific quantita- 
tive measure for the "non-Markovianity" of the dynamical processes. 

We will begin in the following section by a brief summary for our theoret- 
ical framework. On the basis of Kraus representation |15j . we shall explain 
how the dynamics of two independent qubits can be constructed based on the 
single-qubit dynamics. Then in Sees. |ITL| and [TvTI we will apply this formulation 
to examine the non-Markovian dynamics of two qubits subject to, respectively, 
amplitude damping and phase damping noises. Finally, Sec. [V] provides a sum- 
mary of our results and some further discussions. 

II. Theoretical formulation 

Let us consider two spatially separated qubits that are interacting with their 
local environments independently from each other. We shall model the environ- 
ment as a wide spectrum of harmonic oscillator modes so that for each qubit 
the total Hamiltonian reads (we shall take h — 1 throughout) 
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Here wo is the energy separation between the qubit levels |±) (which have ener- 
gies respectively), a z is the third Pauli matrix, and bk is the annihilation 
operator for the oscillator mode with frequency ujk in the environment degrees 
of freedom. The last term Hj in denotes the interaction Hamiltonian for 
the qubit and its local environment, whose explicit form will be specified later. 
Since the two qubits are not interacting with each other and their environment 
noises are uncorrelated, the two-qubit dynamics can be constructed from the 
single-qubit dynamics based on the Kraus formulation, as we shall now explain 

For the single-qubit Hamiltonian 0, suppose one can find the time evolution 
for the reduced density matrix p, it is then possible to express |15j 
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p(t)=J2 E ^ E l ( 2 ) 

i=l 

where p(0) is the initial density matrix for the qubit and Ei are the operation 
elements associated with the time evolution of the qubit. These operation ele- 
ments are encoded with effects of the environment noises on the qubit and satisfy 
Yli=i^l^i — 1 US- Once the single-qubit operation elements are known, the 
time evolution of two independent qubits A, B can then be obtained from the 
operation elements [5] 

Kx = Ef® Ef , 
K 2 = E£® E^ , 
K 3 =E$® Ef , 

K 4 = E$ ® E? , (3) 

where the superscripts A, B are the qubit labels. The reduced density matrix 
p for the two qubits thus evolves according to 
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p{t) = Y J K lP {Q)Kl (4) 

i=l 

with p(0) the initial density matrix for the two qubits. Applying this formula- 
tion, in the following sections we shall consider specific models for the system- 
environment interaction Hj for which the operation elements can be found ex- 
plicitly. We will consider two types of environment noises separately: the ampli- 
tude damping noise, which can cause energy relaxations in the qubits, and the 
phase damping noise, which can lead to dephasing in the qubits PQ. Based on 
these results, we shall study the time evolution of the two-qubit entanglement. 

III. Amplitude damping 

A. Single-qubit dynamics 

As our first case for the qubit-environment coupling, let us consider a single qubit 
interacting with the environment by absorbing and emitting energy quanta from 
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the oscillator modes, and at the same time making transitions between the upper 
and lower qubit levels. This interaction would thus change the level occupations 
in the qubit and cause energy relaxations (thus frequently termed the amplitude 
damping channel). We model this coupling by using the following interaction 
Hamiltonian in (JlJ [TB] 



Hi = (9k<r-bl + g* k <T+bk) , (5) 

k 



where are the coupling constants and a± are the raising/lowering operators 
for the qubit levels. The time evolution of the qubit in the presence of the 
coupling ([5]) can be solved exactly [17] and the corresponding operation elements 
in the basis {|+), |— )} are 



where q = \/T— \p\ 2 . Here p is the solution for the equation 

d p(t) = - f drf(t-T)p(r) (7) 







dt 

with / the noise correlation function 
f(t-r) 

,gl9k'(b k b{,) e <(«o-«-»)t e -iK-« v )T . (8) 

The expectation values in ([8| are taken with respect to the state of the envi- 
ronment. For our calculation, we shall consider a vacuum initial environment 
state. Thus the noise correlation function ^ reduces to 

f(t- T ) = J2\^\ 2 ei{Uo ~" knt ~ T) 

dwJ(w)e l(u »- u)(t - r) (9) 



with J(ui) the spectral function for the coupling |16| . In the following, we 
will consider a specific form for the spectral function that shall enable analytic 
solution for Q. The single qubit dynamics can thus be attained explicitly. 

Let us consider a Lorentzian spectral function for the qubit-environmcnt 
coupling [T3] 

J{u =^7 — a^To To- (10) 

v ' 2vr (uj - cj + A 2 + A 2 v ' 
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Here 7 is the (phenomenological) decay rate for the upper qubit level |+), A is the 
bandwidth of the coupling (thus A -1 characterizes the environment coherence 
time), and A is the detuning from the resonance frequency uiq- Utilizing ( |10[ ) 
in ([9]), one can obtain the correlation function / easily and then solve for p(t) 
from 0. One finds [9] 



, , dt\ X — iA ( dt 
cosh[- + — — smh - 



(11) 



p(t) = e-(^)« 
where 

d = v/(A-iA) 2 -2 7 A. (12) 

With the explicit form for the function p, one can construct the operation el- 
ements in Q. Following the prescriptions of Sec. II. we are now ready for 
studying the two-qubit dynamics. 

B. Two-qubit dynamics 

Let us now consider two non-interacting qubits that are coupled to independent 
amplitude damping noises separately, each described by d5j). With the results 
for single qubit dynamics in the preceding subsection, as described in Se c|II.| one 
can obtain the two-qubit operation elements Ki using ((6} and ^TTJ in ((SfTThe 



time evolution of the reduced density matrix for the two qubits then follows 

from Q. Using the basis {| + +)> | H — ), | — h), | )}, we find the matrix 

elements [9] 

Pn(t) = \paPb\ 2 Pii(0) , 

P2i{t) = bA| 2 (p 22 (0)+g|pii(0)) , 

P33(*) - IPb| 2 0°33(O)+<Za<°1i(O)) , 
Pii{t) = 1 - (pil(t) + P22(t) + P33{t)) , 
Pu(t) = \PA\ 2 PB Pl2(0) , 

Pi3(t) = Pa\pb\ 2 /Ois(0) , 
pu(t) = paPbPu(0), 

P23(t) = PAP* B P23(0), 

P2i{t) = pa (p2i(0) + q 2 B pi3(0)) , 

Pu(t) = p B {p3i(0) + QAPl2{0)) ■ (13) 

Here p a (with a = A, B the qubit labels) are given by ( |11| but with the param- 
eters now carrying the qubit labels, ie. 7 — ¥ 7^, A — > A Q , A — > A a (hence also 
d — >• d a ), and as in (JgJ) , we denote q a = y/l — |p Q | 2 . 

In order to study the entanglement dynamics of the pair of qubits, we shall 
adopt concurrence as the entanglement measure |18j . namely 

C(t) = max I 0, sfjli - yfjii - sfjii - y/pij , (14) 
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where fa (i = 1 ~ 4) are eigenvalues of the matrix p{(jy ® ay)p*(<jy 
in descending orders; Oy is the second Pauli matrix for qubit a. Using (|13| 
and (11), one can find the concurrence evolution for any initial density matrix 
accordingly. For instance, for an initial density matrix of the "X-form" [19] 



P(0) 



/ Pn 







P22 
P32 





pu 

P23 

P33 

p 4 4 



(15) 



it follows from ( 13 ) and ( 14 ) that the concurrence for the reduced density matrix 
p(t) is given by 



C(t) = 2 max {0, \p u (t)\ - VP2 2 (f)p 33 (t) , 

\P23 



(16) 



In the following, we shall present results based on the calculations outlined 
above. We will consider symmetric configurations in which the two qubits are 
identical, so that they have the same level separation uj and the same decay rate 
7.4 = 7b = 7, and are subject to amplitude noises of identical characteristics, 
thus Xa = As = A, Aa — As = A, and consequently pa — Pb — P in (13). 

Fig. [I] displays our results for the entanglement evolution of two qubits with 
the initial density matrix 



p(0) = 





(\ 











\ 


1 





1 


1 







3 





1 


1 

















2 
3 


) 



(17) 



which is an entangled mixed state. In Fig. [TJa) we plot the concurrence evo- 
lution at fixed detuning A = O.OI7 for different values of A. It is seen clearly 
that the crossover from oscillatory behavior to monotonic decay occurs when 
the coupling bandwidth changes from small to large values. In Fig. [ijb) the 
concurrence evolution is plotted at fixed coupling bandwidth A = O.OI7 for dif- 
ferent values of detuning A. One can observe that the entanglement evolution 
changes from damped oscillatory behavior to barely damped one as the detuning 
varies from small to large values. Since entanglement preservation is a key issue 
in quantum information sciences, it is of great interest to analyze in detail how 
such transitions arise. Let us first attempt with simple, qualitative reasoning. 

For the results in Fig.[]Ja), when the coupling bandwidth is small (A C 7), 
it corresponds to the limit with each qubit coupled to almost just one single 
oscillator mode. The entanglement dynamics is thus non-dissipative since the 
entanglement can transfer back and forth between the qubits and the oscillator 
modes involved [2UJ [21] • In other words, in the limit of long noise coherence time 
(ie. large A -1 ) the environment memory can reconstruct the qubit entanglement, 
so that the concurrence displays an oscillation that is hardly damped. In the 
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Figure 1: Entanglement evolution of two^cjubits with the initial density matrix (17 1 
subject to amplitude damping noises. We plot the concurrence for (a) different A at 
fixed detuning A = O.OI7 and (b) different A at fixed coupling bandwidth A = 0.0f7. 
In (c) the solid curve depicts the entanglement evolution for large coupling bandwidth 
A = 10. 07 (with A = O.OI7), and the dashed curve for large detuning A = 10. 07 (with 
A = O.OI7); note that for clarity, the time axis in (c) has a smaller range than those 
in (a), (b). The axes A in (a) and A in (b) are in units of 7, while the time axes in all 
plots are in units of I/7. 



other limit of large coupling bandwidth (A > 7; see the solid line in Fig.Qc)), the 
qubits are coupled to a wide spectrum of oscillator modes. Thus the environment 
memory time becomes exceedingly short, so that the entanglement dynamics 
becomes Markovian and the concurrence decays to zero monotonically. In the 
intermediate regime, due to the finite environment memory time, concurrence 
oscillation is damped but approaches zero in a non-monotonic manner. Since 
entanglement oscillation is an indication for non-Markovian dynamics [8lfTHl9"]. 
the crossover associated with the coupling bandwidth can also be regarded as 
one between the limits of Markovian and "extremely" non-Markovian dynamics. 

For the results in Fig. [TJb), we note that for large detuning (A > 7; see also 
the dashed line in Fig. [TJc)) the qubits can only couple weakly to the environ- 
ment noises. The concurrence is thus barely damped in this limit. With de- 
creasing detuning, the environment noise can couple more and more effectively 
to the qubits, which results in stronger and stronger damping in the concur- 
rence evolution. Therefore, the crossover from undamped to damped oscillatory 
entanglement dynamics that we are seeing in Fig. [ljb) is connected with the 
transition between the weak and strong coupling limits in the qubit-environment 
interaction. 

In order to substantiate the arguments above, let us now turn to more quan- 



titative, detailed analysis for the behavior of p(t) based on ( 11 ) 



Small coupling bandwidth 



In the limit of small A, suppose A <C A <C 7, one can find from ( 12 1 that 

d ~ iy/2-yX + A 2 =iuj d . 



Hence it follows from (111 
p(t) ~ 



cos 



(18) 



(19) 



which oscillates in time with no damping. This gives rise to the undamped 
concurrence oscillations, which are manifestations of non-dissipative dynam- 
ics. Note that this regime is small in Fig. [IJa); it can be seen more clearly in 
Fig. [ljb). For larger bandwidths A<A<7, one can obtain from ( 12 1 



and thus (11) would yield 
P(t) 



e 2 cos 




(20) 



(21) 



We see that p(t) is now oscillating with frequency \J 7A/2 with a damping time 
~ A -1 . Therefore for A <C 7, there can be many oscillations within the damping 
time; one would thus find weakly damped concurrence oscillations in this regime. 
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It should be noted that in the limit under consideration, if the detuning is 
large, such that A 3> 7 3> A, then we have from ( 12 1 



d~\-iA. (22) 
It then follows readily from ( |TT| ) 

p(t) ~ 1 . (23) 

The concurrence thus remains a constant throughout the time evolution; see the 
dashed curve in Fig. [ljc). As explained above, this is because the detuning is 
so large here that the qubit-environment coupling becomes no longer effective. 
The qubit dynamics thus becomes entirely immune to the environment noise. 

2. Large coupling bandwidth 

In the limit of large coupling bandwidth, suppose A 3> 7 3> A, it follows from 



([12j) that 

d~A-7-iA. (24) 
One can thus obtain from (JTTJ) 

p{t)~e-i\ (25) 

which is typical Markovian behavior. The concurrence evolution is therefore 
strongly damped in this regime and decays to zero monotonically (see the solid 
curve in Fig. |ljc)) . Note that in arriving at (25) the detuning has been taken 



to be small, ie. A <C A, 7. If we have large detuning, so that A ^> A 3> 7 or 
A 3> A 3> 7, it is easy to show that one would get p(t) ~ 1 in both cases. The 
concurrence thus becomes stationary due to the ineffective coupling at large 
dctunings. 



3. Spectral function 



To further justify the foregoing analysis in 1. and 2., it is also instructive to 
examine the spectral function J(ui) in the corresponding limits. In the limit of 
small coupling bandwidth, one has from ( 10 ) 



7A 

lim J(oj) = — S(oj — ujq + A) 

A->0 2 



(26) 



which signifies that, as noted above, each qubit is coupled to a single oscillator 
mode with frequency ujq — A. As a consequence, a non-dissipative, oscillatory 
concurrence evolution follows. In fact, making use of (26 1 in j9j), and then 
solving p(t) from Q, one can recover (19) accordingly. In the limit of large 
coupling bandwidth, one finds similarly from ( 10 1 



lim J(oS) = — 

A— >oo 27T 



(27) 
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Therefore, we have a white noise spectrum for the qubit-environment coupling 
in this limit. The entanglement dynamics thus becomes Markovian and the 
concurrence decreases to zero monotonically. Again, using (27) in j9j, one 
can regain (25 1 easily from Q. In the intermediate regime between the two 
extremes, the concurrence displays damped oscillations, where the oscillations 
are manifestations of the finite environment memory times due to non-zero 
coupling bandwidths. 

As another example, we consider an entangled pure initial state which has 
the density matrix 



P(0) 
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(28) 



The results for the calculation are shown in Fig. [2] Although there are now 
regions with entanglement sudden death (and rebirth), our explanations above 
for the concurrence oscillation/damping in Fig. [I] can be carried over. 



IV. Phase damping 
A. Single-qubit dynamics 

Let us now turn to another class of qubit-environment interaction that couples 
longitudinally to the (pseudo)spin of the qubit. The environment noise would 
therefore not cause any energy relaxation in the qubit, but randomize the rel- 
ative phase between the qubit levels. In the single-qubit Hamiltonian ([lj, this 
coupling can be modeled with the following qubit-environment interaction |22j 

ff/=$>*(fl4 + - (29) 



where, as previously, gu are the coupling constants. Unlike ([5]), here Hj com- 
mutes with the bare qubit Hamiltonian (u;o/2)cr z , and thus would not induce 
any population transfer between the qubit levels |+) and |— ). The coupling 



(29) would instead act like a fluctuating magnetic field along ^-direction that 
can randomize the phase of the qubit 22J. This qubit-environment coupling is 
thus often referred to as the phase damping channel pQ. 



For the interaction Hamiltonian (29), one can work out the time evolution of 



the qubit explicitly [22] . It is then not difficult to find the corresponding single- 
qubit operation elements. Suppose the environment is initially in the vacuum 
state, in the usual basis {|+), |— )} one can find 

Ms !)■ Ms »)• (30) 

where, as before, q = \J\ — \p\ 2 . Here we have 

p(t) = exp{r(t)} (31) 
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Figure 2: Same as Fig. [T] but for the initial density matrix (28 1. The flat regions at 
the base planes of (a) and (b) are signatures of entanglement sudden death. 



with 



r(i) = -2^| 5fc | 2 f dr [ dT\ 

Jo Jo 



—iUkir-T 1 ) 



(32) 



a real negative function of time. As in the case of amplitude damping, we 



can rewrite T(t) of (32) in terms of the noise correlation function f(t) and the 
spectral function J{lo) using we get 



r(t) = -2 f dr f dT'f(T-T')e- iuo{ - T - T "> 
Jo Jo 

duj J{uj)e- lu{T - T ' ] 



-2 f dT I dr' 
Jo Jo 



(33) 



As previously, once a specific form for J(uj) is given, it is then possible to find 
the qubit dynamics accordingly. 



Let us consider again the Lorentzian spectral function (10 1. Note, however 



that since the interaction ( 29 1 would not cause any transition between the qubit 



levels |+) and |— ), the frequency ojq no longer plays the role of the "resonant 
frequency" . One can thus redefine lj c = uq — A as the "central frequency" for 



the coupling (this follows from the form of the spectral function (10), which 
centers at uj = oj c ). Applying (10) in (33), one can obtain 



r(t) = -2 7 a< 



A 



A 2 + 



1 



(A 2 + w 2 ) 2 



(A 2 - wjj) (1 - e~ A * cos(w c i)) 



+2Au; c e A * sin(w c t) 



(34) 



It is interesting to note that for non-zero ui c , T(t) (and hence p(t)) can oscillate 
in time. The qubit decoherence can consequently occur non-monotonically. In 
the case of two qubits, this can lead to entanglement oscillations, as we shall 
now show. 



B. Two-qubit dynamics 

For two independent qubits A, B subject to uncorrelated local phase noises, 



again, we resort to the prescription outlined in Sec. II. Using (30) in ^ 



one can derive the two-qubit operation elements in the phase damping channel. 
Substituting the results into Q, one can then obtain the time evolution of 
the the matrix elements for the two-qubit reduced density matrix (in the basis 
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->}) 




Pii{t) = 


p H (0), z = l~4 


Pu(t) = 


Pb Pu(0) , 


Pis(*) - 


PA (013(0) , 


= 


PAPB Pu(0) , 


023 (*) = 


PAPB P23(0) , 


P2i{t) = 


PA P2A(0) , 


P34{t) = 


PB P34(0) , 



(35) 

where p a are the expressions for qubit a = A, B in accordance with ( |31[ ) and 
(34) (with appropriate qubit labels added). Note that here we have made use 



of the fact that pa,b{^) are real functions of time. From (35 1 we see that, as 
noted previously the phase noise does not affect the qubit occupations since the 
diagonal elements of the reduced density matrix remain unchanged. However, 
it diminishes the coherence of the density matrix by corroding the off-diagonal 
elements of the reduced density matrix. 



In order to examine the entanglement dynamics ensued from ( 35 1 , again 



we shall consider initial density matrices of the X-form (15). The concurrence 



of the qubits can then be calculated using (31), (34) together with (35) in 



(16). As previously, we shall consider symmetric configurations where the two 
qubits are identical (so that they have the same level separation ujo and the 
same decay rate 7a = 7s = 7) and the phase noises acting on them have the 
same characteristics, namely Xa = Xb = A, ujf — u)f — uj c , and consequently 
Pa =Pb =P in d35b . 



Fig. [3] illustrates our results for the entanglement evolution of two qubits 
with the initial density matrix (17). In Fig. [3]ja) we plot the concurrence evo- 
lution at fixed central frequency uj c = I.O7 for different coupling bandwidths 
A. When the coupling bandwidth varies from small to large values, crossover 
from oscillatory behavior to monotonic decay in the concurrence evolution can 
be seen clearly. In Fig. [3jb) the concurrence evolution is plotted at fixed cou- 
pling bandwidth A = O.OI7 for different central frequencies uj c . One can observe 
that the entanglement evolution changes from monotonically decaying behavior 
to barely damped oscillations as u> c varies from small to large values. Another 
set of results for the entanglement dynamics under phase damping noise is il- 
lustrated in Fig. [4j where we consider the initial density matrix ( 28 1 . Again, 
similar crossovers between different regimes are clearly visible. 

Qualitatively, the results in Figs. [3] and [4] can be understood in the same 
manner as presented in Sec. 



III. 



for the entanglement evolution under ampli- 
tude damping noise. For instance, in both Figs. [3ja) and |4][a) the regime with 
small coupling bandwidths correspond to the limit with each qubit coupled 
to one single oscillator mode. The entanglement dynamics thus exhibits non- 
dissipative behavior, as entanglement swapping can take place between the pair 
of qubits and the two oscillator modes. The undamped entanglement oscilla- 
tions in this regime are consequences of the long environment memory time A -1 , 
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Figure 3: Entanglement evolution of two^aubits with the initial density matrix (17 1 
subject to phase damping noises. Here we plot the concurrence for (a) different A at 
fixed central frequency cj c = I.O7 and (b) different ui c at fixed coupling bandwidth 
A = 0.0f7. In (c) the solid curve illustrates the entanglement evolution for large 
coupling bandwidth A = 10. 07 (with uj c = I.O7), and the dashed curve for large 
central frequency ui c = 10. 07 (with A = O.OI7). For clarity, we plot in (c) a smaller 
range of time than in (a), (b). The units for the axes are the same as in Fig. [l] The 
flat areas at the base planes in (a) and (b) indicate regions with entanglement sudden 
death. 



which enables entanglement reconstructions in the qubits. For large coupling 
bandwidths, the wide spectrum of the qubit-environment coupling entails an 
extremely short environment memory time. Therefore, the Markovian entan- 
glement dynamics emerges and the concurrence decays to zero monotonically 
(see the solid curves in Figs.[3jc) and[4]jc)). For the intermediate regime, the fi- 
nite coupling bandwidth implies non- vanishing environment memory time. The 
concurrence thus carries out damped oscillations, in which the oscillatory part 
reflects the memory effect that survives the phase damping noise. 

As pointed out earlier in the present section, the central frequency w c is 
the peak frequency for the qubit-environment coupling. Therefore, at large w c 
the qubits are coupled to the phase noise only very weakly. This is why in 
both Figs.gb) andgb) the concurrence is barely damped at large u> c (see also 
the dashed curves in Figs. [3^c) andQc)). With decreasing uj c , as the noise 
can couple more and more effectively to the qubits, damping in the concurrence 
evolution gradually sets in. Eventually, in the limit of small w c , the entanglement 
oscillation disappears and the concurrence drops to zero monotonically. 



As in Sec. III. more quantitative analysis for the results here can be achieved 



by looking into the behavior of p(t) according to (31 1 and (34 1 in the limits of 
small and large coupling bandwidths. It should be noted, however, that unlike 
amplitude damping noise, here the regimes of small and large coupling band- 
widths are determined from the relative magnitude between A and uj c (instead 
of the qubit decay rate 7). This is because under the action of phase damping 
noise, the qubit populations would not change. Thus the time scale 7 becomes 
inessential to the non-Markovian dynamics of the qubit (see later). Instead, it 
is the central frequency lo c that comes into play (which is clear from the explicit 
form of T(t) in (34)). Therefore, for phase damping noise although the central 
frequency uj c plays a role similar to the detuning A in amplitude damping noise, 
the concurrence oscillation depends more sensitively on w c (compare Figs. |T|(b) 
with[3]jb), and Figs. [2jb) with Eft)), especially for regions with small A and 
w c ). With the above observations, let us now examine p(t) for small and large 
coupling bandwidths using (31) and (34). 



1. Small coupling bandwidth 



For narrow coupling bandwidths A -C oj c , we find from (34) that 



T(t) 



2 7 A 



[l + At - e~ xt cos{u c t)] 



(36) 



Therefore for times t <C A -1 , pit) = exp{r(t)} oscillates with negligible damping 
and the concurrence dynamics exhibits non-dissipative characteristics. This can 
again be anticipated from the spectral function J(uj), which reduces to the form 
for single mode coupling in this limit (see (26)). Indeed, using the spectral 
function (26) in (33), one can derive an expression for T(t) which has exactly 
the same form as (36) except with Xt = (ie. being totally non-dissipative). For 
Xt ~ 1, however, we see from (36) that damping begins to set in. As a result, the 



oscillating cosine term in (36), which represents the environment memory effect, 
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is attenuated and a damped concurrence oscillation thus follows. At long times 
t ^> A -1 , one can infer from (36 1 that p(t) eventually drops to zero exponentially, 
thus also the concurrence (save the complications due to entanglement sudden 
death, such as in Fig. |3| . 

Note that in the present limit, if the central frequency is so large that u> c 3> 



7, A, it follows from (36) that T(t) would become vanishing. In this case we have 
pit) ~ 1 and the entanglement becomes stationary, as the qubit is completely 
decoupled from the environment in this limit. This can be seen from the dashed 
curves in Figs. [3^c) and[I]jc). 

2. Large coupling bandwidth 



For large coupling bandwidths A 3> u> c , it follows from (34) that 

~ xt cos(w c t) * 



r(t) 



-2 7 t- 



1 - e 



X 



(37) 



Since A is now large, the damping time (or environment memory time) A -1 
is short. Within the damping time, the oscillatory cosine term in (37) cannot 
even complete one period since w c A _1 <C 1 here; the environment memory is 
therefore washed away quickly. For times t A -1 , we get from (37) that 



P (t) = exp{r(t)} 



-2ft 



(38) 



which manifests the Markovian nature of this limit. As a consequence, the 
concurrence becomes strongly damped and exhibits Markovian behavior [TTj . 
As previously, we can further corroborate the argument above by considering 
the spectral function J(u>) in the present limit; the white noise spectrum (see 



Eq. (27)) is again responsible for the Markovian entanglement dynamics in this 
limit. As a check, one can use (27) in (33), which would yield (38) immediately. 



V. Discussions and Conclusions 

In summary, based on exactly solvable models, we have studied the entangle- 
ment dynamics of a pair of non-interacting qubits that are separately coupled to 
their local environments. In particular, we examine in detail the crossover be- 
tween the limits of dissipative and non-dissipative entanglement dynamics that 
are manifested in the damped/undamped concurrence oscillations. We show 
that this crossover is connected with the bandwidth of the qubit-environment 
coupling. We have also studied the crossover between the limits of strong and 
weak couplings, which is controlled by the detuning of the qubit-environment 
interaction. Both categories of crossovers have been considered for two types 
of environment noises, the amplitude damping channel and the phase damping 
channel, and physical pictures for the results have been provided. 

Although we have illustrated our results only for two specific initial density 



matrices (17) and (28), our conclusions are indeed fairly general. For instance, 
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suppose only one of the qubits is coupled to the environment, it is easy to derive 
from our formulas that the concurrence evolution for any pure initial state would 
follow the factorization law [23] [9] 

C(t) = \p(t)\C(0), (39) 



where p(t) is given by (11) and pl| ), respectively, for amplitude and phase 
damping noises. Since for maximally entangled pure initial states C(0) = 1, 
it follows that p(t) completely determines the concurrence evolution of the the 



maximally entangled state for this noisy channel. Therefore, according to (39), 
for any pure initial state the concurrence evolves according to the action of the 
noisy channel over the maximally entangled state [23] ■ The oscillation/damping 
of p(t) in different regimes would thus exhibit also in the concurrence evolution 
for general pure initial states. Similar results exist also for certain class of mixed 
initial states when both qubits are subject to local noises [9]. 

An issue that we did not address in this paper is the interplay between am- 
plitude noise and phase noise. For Markovian entanglement dynamics, Yu and 
Eberly had demonstrated that the combination of phase damping and amplitude 
damping noises can give rise to non-additive entanglement dissipation [24] , For 
the non-Markovian case, it is certainly of interest to investigate whether addi- 
tional effects would emerge when both phase and amplitude noises are present. 
Also, in our analysis of the entanglement evolution, we did not examine the de- 
tailed dynamics close to entanglement sudden death. We wish to take up these 
issues in future works. 
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